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Ai^STKACT 



Pile oscillators have been successfully emi^loyed to detect 
small concentrations of in;purities in substances such as boron in 
graphite. In several instances the technique has proved more 
sensitive than the commonly used "danger coefficient" method. The 
danger coefficient method using the Aerojet General Nucleonics 201 
Reactor has been demonstrated to be as sensitive as the pile 
oscillator method on other reactors; thus a pile oscillator for the 
AGN 201 Reactor offered the possibility of a further improvement. 
Because of low cost and simplicity, a design utilizing one of 
the installed ion chambers in the reactor and the linear amplifier 
on the control console was chosen for construction of the pile 
oscillator. Additional instrunientation employing the output of the 
linear amplifier gave excellent performance in reproducing a test 
signal; however, when the pile oscillator was put into operation 
satisfactory reproduction of the output for a given sample could not 
be obtained. It is felt that the error is due to reactor noise. 

It is concluded that a pile oscillator utilizing the installed 
ion chai.ber and linear amplifier of the AGN 2C1 Reactor operating at 
its present power level is not satisfactory because of reactor noise. 
It is doubtful that any substantial improvement could be obtained 
unless the operating power were increased above the presently 



licensed value. 
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!• Introduction. 



The time behavior of a nuclear reactor i. ay be affected to 
varying degrees by the introduction of a material into the reactor* 
This behavior can be related to the nuclear cross sections^ of the 
material. Methods have been devised to utilize the time behavior 
characteristics of the reactor to measure, in particular, neutron 
absorption cross sections. 

Certain important impurities have microscopic absorption cross 
sections which are very large. Siaall percentages in materials will 
increase the total macroscopic cross section significantly. Boron is 
a common exoi.ple and its presence in graphite or iron must often be 
quantitatively deteniiined. The methods by which absorption cross 
section measurements are made may readily be adapted for measuring 
the concentration of impurities. This is usually done by a comparison 
with a standard. 

One technique is called the danger coefficient niethod 1 | ^ . A 
sample of material is placed into the reactor in a fixed position, 
usually the center of the core. The time behavior of the reactor is 
noted and through appropriate relationships a comparison made with a 
standard sample. 

Another method, first used by Langsdorf pil® 

oscillator. A mechanical drive oscillates a seunple within the reactor. 

I 

A neutron-sensitive ion chamber is used to record the overall power 
of the reactor. The oscillating seunple causes an alternating 

^Elementary definitions of the terms and concepts which are taken 
fro:, nuclear physics and reactor theory are given in Appendix I. 

2 

Numbers in brackets refer to the Bibliography. 
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fluctuation in the overall reactor flux. The overall flux containing 
the variational couiponent is sensed by the ion chainber. This signal 
current is treated so that the variational component is rectified and 
integrated with respect to time. The output voltage of the integrator 
is calibrated with a standard seycple permitting the use of the device 
for measuring either cross sections or concentrations of impurities. 

In addition to the variation of overall flux caused by the 
oscillating absorber, there is a significant local variation. By 
placing the ion chamber close to the absorber the local signal can 
also be detected. If thenral neutron cross sections are being 
measured, the best location for the oscillator and ion chamber is in 
a region such as the graphite reflector where neutrons are predominantly 
themial. In this manner the overall reactivity changes induced by the 
sample are insignificant compared to the local disturbance. This 
technique was described by Hoover, et. al. ^3^ . The sample was 
oscillated through an annular ion chamber located in the graphite 
reflector of the X-10 reactor at Oak Ridge. The output of the ion 
chamber was aii:plified by a special feedback amplifier which passed 
only the variational component of the overall signal. This output 
was fed to an integrator after appropriate rectification. This method 
has been successfully used at Batelle Memorial Institute ^4^ • 

Comparing the pile oscillator method to the danger coefficient 
method, Langsdorf reported an increase in sensitivity of a factor 
of 2C over danger coefficient measurements with the same reactor. 

Using the Aerojet General Nucleonics 201 Reactor (AGN 201) for 
danger coefficient measurements, Carew and Mow made measurements 

as precisely as had been done with other reactors using pile 
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oscillators* It therefore appeared possible that a pile oscillator 
for the AGN 201 might produce even more precise measurements. 

It was decided to construct a pile oscillator for the AGN 201 
reactor and investigate its characteristics* Many possible designs 
had to be considered. In oriler to oscillate a sample through the 
core of the reactor, the glory hole iv'ould have to be used. The 
physical construction of the reactor would preclude plo,cing a 
detector close enough to detect the local flux variation. If the 
oscillator were located in a region where a close-in detector could 
be placed, the access ports or thermal column would presumably be 
used* In this case the overall reactor flux would not be greatly 
affected and the reactor would be acting primarily as a flux fi^ld. 
All the previous devices of the type that oscillated within a thermal 
column or reflector were used in reactors having a thermal flux 

9 _2 -1 

of the order of 10 neutrons cm sec * Since the thermal flux of 

6 - 2-1 

the AGN 201 is of the order of 10 neutrons cm sec , this method 
would be at a disadvantage* 

Thus the type of oscillator which measured the changes in the 
overall reactor flux was indicated* From considerations of simplicity 
and reasonable expense, a design was base i on the employment of one 
of the penaanent ion chnnibers of the reactor and its associated 
amplifier on the c'ontrol console* It was felt that the following 
procedure would be appropriate: 

1. To construct a mechanical drive mechanism which would 
oscillate a sample of material within the reactor with 
exactly reproducible motion over each cycle and with 
accurate variable speed control . 
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2, To devise an inst riiii entation system i/hich would reproduce 
accurately an input which simulated the actual signal 
resulting frorr the oscillation of a typical sample. 

3. To perform experir.ents oscillating samples v/ithin the 
reactor. Using the instrumentation system devised, to perform 
measurements to determine if t)ie output reading of the system 

, is reproducible and to determine the sensitivity of the 
overall technique . 

2. Reactor kinetics for an oscillating absorber. 

The complete fon.ulation of the effects of an oscilleting 
absorber in a reactor renuires extensive developi ent and involved 
expressions, i rigorous treatment is available from Weinberg ^ j • 
However, since the pile oscillator method to be used with the AGN 201 
reactor measures overall changes in neutron flux, the j)roblem is 
considerably simjilif ied. Since interest is therefore confined to the 
tiiewise variation of flux rather t^san the spacewise variation, the 
case L.ay be considered in which the time-dependent part of the 
solution of the governing equations is separable from the geo?, etry- 
dependent parts* This is achieved if the reactor is at or near the 
critical condition. 

In the following treatment the assui ption of a homogeneous 

reactor is made. The one-group space-independent kinetic equations 
1 

are; 

^See Appendix I for developL.ent of these equations and definitions 
of syi:bols. 



4 



6 



( 1 ) 



dn Ak (AIc + 1) ^ ^ ^ „ 

r'T"- p — Z-. ^i'^i 



dt 



dC 



i = 1 



n- A. C. 

dt '1 11 



( 2 ) 



In order to study the effect of an oscillating absorber, let it 
produce a change in k which is sinusoidal in form. Writing this as 

k(t) -Ak,, (j2 . -1) 

and limiting to small values,^ small amplitude 

variations of n(t) and C.(t) are 



jcjt 



n(t) = n -f n, e 
' ' o 1 



C. (t) = C. + C., e 

i' ' 10 il 



jot 



Mrhere ■ ■ 1 and n is the 

n o 

o 



neutron density in the critical state. Substituting these expressions 
and corresponding derivatives into equation (l) gives 



a> 



n, ej^^=:^n e^^ ^ n. -.^n 0^^*- £ n 

J? (D I Q) i ls> 



P* 2,jcJt A .ioi- 
: — n, e - P n, e^ + 



(3) 



i 1 






T ^ 



a> 




C. + \ ^.C.e 

' 1 10 / i 1 



jcjt 



i = 1 



If terms and are neglected as being second order, and term (5) 
is neglected as being small in comparison vith term (g) , the equation 
reduces to 



6 



H 



jcon^ :;o - H «1 + N Xi 

£ ^ 



( 4 ) 



i = 1 



Repeating the process for equation (2) gives 
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( 5 ) 



t Pi 



. „ jiot ^*^0 jot 2jo 

JCjCii eJ = _ eJ + _ e ^ + J" "<> 



P i jc^>t ^ n ^ r. 

+ z A«C. +A.C,|0 

/I 1 10 1 11 



Again neglecting second order terii a gives 



JcJC., cr 



(5. „ + -^n, - >. C.. 

d ^lO ^llJ 



( 6 ) 



Solving this equation for C^j and making the substitution into 
equation (4) gives 



jconj 



^ o _ 




6 . 


■Pi "i 1 


6 




1 

o 


T"* 


i = 1 


Jio+ X . 


i = 1 


jto+ A. 
*. * 



( 7 ) 



Discarding the last term under the assumption ^ 



solving for ^1 
n 




n 

o 



AK 



joji+ (3 



i 

i = 1 



( 8 ) 



= -y" W(j60) 

where W(j6j) = 




( 9 ) 



Thus it has been shown, subject to the approximations made, 
that tlie ratio of the variational component of the neutron level 
to the average neutron level is proportional to the input change in 
reactivity for a given value of CO* W(jcj) is called the reactor 
transfer function. 



* 
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The results of a set of calculations for different parameters 
, and for varying CJ is given in reference | 7^ . The particular 
results of interest and applicable to the AGN 2C 1 reactor are shown 
in Figure 1. 
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CJ (radians/sec) 

Amplitude of W(jtj) as a Function of cj 
Figure 1 
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3. Description of pile oscillator for the AGN 2C 1 reactor, 
(a) General arrangement. 

Figure 2 shows a block diagrain of the pile oscillator 
constructed for the AGN 2C1 reactor. 




Figure 2 

Block Diagram of Pile Oscillator for AGN 201 

The mechanical drive (l) moves the sample of material (2) 
through the core of the reactor with uniform linear motion. The ioti 
chamber (3) detects the overall neutron flux level including the 
change induced by the saiuple and su])plies a signal current to the 
linear amplifier (4) installed in the reactor control console. The 
output of the linear amplifier goes to the first sto.ges of the 
metering system (5) where the signal representing the response to the 
perturbation of the oscillating sarple is removed f roDj the overall 
signal. Tlie positive part of this voltage is integrated with 
respect to time by the integrator (G) froii which the output may be 
read directly on a voltmeter or recorded on a strip chart recorder (7) 
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for graphical analysis. The output voltage of tiie integrator is the 
criterion for the comparison of an unknown sample with a known sample, 
the frequency and number of cycles being fixed. 

(b) Mechanical components. 

In the design of the mechanical drive it was desired to provide 
an arrangement which would be reasonably simple to manufacture and 
which would be adaptnble to changes found necessary in the course of 
the investigation. It was also desired to drive the sample in such a 
manner that the wave fonr of the induced change in tlie multiplication 
factor would be sinusoidal. Since the absorption ic:portance function 
of the reactor (see Appendix II) can be approximated by a cosine 
squared function, a linear movement of the sample is equivalent to a 
sinusoidal sarple motion through a linear flux gradient. Thus the 
mechanical drive was designed to produce linear movement of the 
sample through the reactor core. 

The drive mechanism is shown in Figure 3. A constant speed 
electric motor drives a hydraulic speed reducer to which is connected 
one sj)rocket for a roller cliain. Another sprocket is mounted so as to 
stretch a heavy roller chain (3/4 in. pitch) between the sprockets. 

The connection of the chain is i.:ade by a tool steel pin with a one 
inch overhang. This pin engages a vertical brass guide slot on an 
aliurinum block whicii is free to move horizontally on two guide rods. 
Thus the travel of the pin causes the block to follow in a linear 
motion except at the end points where the direction of motion is 
reversed. The sample holder is attached to the sliding block. 

The sajTiple holder is made from 3/4 inch o.d. aluninuit tubing 
having a wall thickness of j/32 inch. The tube was filled with 
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Figure 3 
Mechanical Drive 
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reactor grade grapnite which waft nac’iined to close dimi.etral 
tolerances. A portion of the aluL.inui., tube was cut away at the 
iridpoint. An interruption in the graphite at this point permitted 
the introduction of a tost sample. The graphite filler was employed 
to }. inii ize flux variations caused by the scattering of neutrons by 
the s<ur.ple. 

(c ) Instrumentation • 

The instrumentation system converts the variations in neutron 

flux into useful information. The signal is initiated in t!ie ion 

choi'ber in the reactor and amplified by the linear airplifier on the 

reactor control console. Because of the low power level of the reactor 

the ion chamber current at rated recactor power of one hundred 

-9 

milliwatts is only of the order of 1C ’ amperes. The linear amplifier 
is used therefore to increase the power level signal to usable levels. 
It is fron: this point that the additional instrumentation tahes its 
information. In spealiing of the **signal’* the reference will be made 
to the output of the linear amplifier. 

The signal voltage consists of the s\im of three parts: a D.C. 
voltage representing the average reactor power level, an oscillating 
component representing the response to the j)erturbation of the 
oscillating sairple, and noise. Since the oscillating co. ponent is one 
or more orders of magnitude lower than the D.C. component, it is 
necessary to se])xrate it Croii- the overall signal to permit measurement 
The frequency of oscillation is rpiite low u.aking removal of the 
variational component by usual amplifier techniques difficult. 

An an*plifier of the design used with the Oak Ridge oscillator ^ 3 "^ 
was constructed. This perforued well with low f reouencv sinusoidal 
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inputs but was not sufficiently stable when the inpi|t voltage 
contained a significant an ount of noise. This was eventually 
abandoned. 

The problem was satisfactorily solved by eir.ploying operational 
ani|)lifiers of the type employed in analog computer applications. A 
systei. was developed which subtracted the D.C. level f ron: the overall 
signal leaving the variational voltage. This voltage was rectified 
and integrated with respect to time. The integrator voltage after a 
selected nuiiiber of oscillations was the basic value from the 

4 

measurement. 

In order to describe the instrumentation system in detail it 
will be necessary to first describe two elenentary ann.log computer 
c i rcui ts. 

A high gain operational ajuplifier is symbolized by 
If the output is fed back through a resistor to the input and a 
voltage ej is fed to t!ie input through a resistor Rj : 





the output voltage will be 
Rr 



e = - 
o 



R 



1 J 



= - G 



wh^re G = ~ 



1 



If adiitional voltages are fed into the input in like Lanner, the 



negative sum 



e = - 
o 



R. 



R. 



Rj ^ 



will be formed. 




If U- = R., = etc, G is constant; the shortened symbol 

I 2 d 

will be used for this case. 

If a capacitor of capacity C is used for feedback in place of 
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the resistor, 




e , t!»e integral e =- ^ e, dt 



will be formed. Likewise, additional inputs may be included 



and 




dt where G . = — — 
1 R.C 
1 



If all the input resistances are enual , the shortened symbol — 
will be used* 

It was stated previously that the D.C. component of the signal 
voltage was to be found in order to subtract it from the overall 
voltage leaving only the variational component. This requires a 
definition of the variational voltage. This may be done in the 
following manner: 

I e dt ^ e 
I var var max 

That is, the time integral of the variational voltage will not grow 
with time but will remain within the limit given. This may be 
expressed functionally: 

,t 

e 





- 1 I (e - ej^) dt £ 

RC 



var 



RC 



max 



where e is the total signal , and e^ is the instantaneous value of the 
D.C. component. In other words, if the integral of (e - e^^) is formed, 
the value of the integral will not grow with time but will remain 
within the limit given. 

With this idea in mind the following scheme suggests itself 




e cr -e, 
o b 
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Although the behavior of this system can be described in precise 



terminology, it ray be understood by considering the effect if e is 

o 

not equal to the defined e^« If is nurrerically too small, there 
will be a net positive voltage input. This will generate a larger e 
but, since this is negative by the sign reversing characteristic of 
the injt/egrator , the net positive input will be din.inished until 
equilibriui; is established. 

Thus far no account has been taken of the possibility that the 
average power level of the reactor may be drifting. This would have 
the effect of making e^ drift. This can be minimized by selecting 
proper values of resistance and capacitance in the integrating 
oTuplifier so that it will follow the change in the average voltage. 
Generally speaking, the shorter tiie time constant of the integrating 
oiplifier, the better it will follow the drift. However, if the time 
constfUit is too short, it will tend to follow the oscillating 
component of the input as if it were a drifting D.C. Thus a 
compromise is reached based on the maxir.mn drift anticipated. 

With the D.C. component determined, the variational component 



can be found: 




In order to render the effect of drift insignificant, two more 

amplifiers was added repeating the above loop. The variational 

voltage is then rectified by a diode leaving only the positive halves 

which is the information signal, e. ^ . The use of half wave 

rectification insures that the part of the signal which becomes 

e. ^ corresponds to the sample position being in that part of the 
info o I 
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reactor vrhere the cosine waveforn is closely approximated (see 
Appendix II). The rectified voltage, input to an 

integrator. 




In order to check the reproducability of this system, a 
simulated input was formed. This consisted of an electronic function 
generator with a battery in series* This gave a D*C* voltage with a 
sinusoidal variation superimposed. The magnitudes and frequencies 
were chosen to correspond to the output of the linear amplifier. The 
system was allowed to reach equilibrium such that it balanced the D.C 
cor?ponent and then integrations were performed for a fixed number of 

I 

cycles. This process was repeated, and the resulting output voltages 
of the integrator were read. The readings were always the snme within 
a fraction of a percent. A slightly drifting voltage was added with 
no noticible change. 

A schematic diagram of the amplifier system is given in Figure 4 
Illustrations of the physical arrangement are shown in Figure 5. 

A strip chart recorder was included in the instruji^entation so 
that the final integrator output or any intermediary voltage could be 
recorded. 
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Figure 5 

Instrumentation Equipment 
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4, Operation of the ])ile oscillator. 

In order to Uia\e a measurement, a standard specimen of iron was 
inserted in the sai:ple rod and the rod positioned in the drive 
mechanism so that the sample was in the center of the reactor core 
when the drive mechanism was in its i id-position. The amplifiers in 
the instrujn.ent system were balanced and all inputs and condensers 
shorted. The drive mechanism wa-i started and the reo.ctor control rods 
.idjusted until the trace on the linear recorder on the control 
console showed no tendency to drift. This required several minutes# 

A.11 shorts were removed from the instrumentation system except for the 
final integrating capacitor. AXter several minutes the instruments had 
steadied out and an integration could be coimenced. This was done when 
the drive mechanism was in an extreme ”out” position and continued for 
a fixed number of cycles after which the output of the integrator was 
recorded. The integration for the some sample was repeated for 
approximately five readings. Then the seui^ple was removed and replaced 
by a sample having an unlcnown percentage of impurity and the procedure 
repeated. In this manner the unicnown sar.iple could be related to the 
standard sample which was for all purposes free of the impurity. 

Several difficulties are inherent in the system which limit the 
total time allowable for the measurement of each sample. It is 
difficult to maintain the average level of the reactor at precise 
criticality. The nature of the instrumentation system is such that it 
must be balanced at least every half hour. The choice of component 
values in the instrumentation was to a large part governed by the 
ability to compensate for drift effects. This in turn required that 
the oscillating frequency should not be too low or the instrir ents 
AV^ould tend to follow t ’ie oscillation as if it were D.C. drift rather 
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than the information bearing signal. 

The choice of oscillating frequency must take other factors 
into account. From the plot in Figure 1 it would appear that a lower 
frequency of oscillation is desirable. However, statistical 
considerations raise the possibility of advantage from higher 
frequencies since a greater number of observations can be acciir.iulated 
in the same total tiii.e. This argnu ent assai.es that the relative error 
for each cycle does not increase as the frequency of oscillation is 
increased. Weighing all the contrasting factors it appears that a 
frequency of oscillation near the inflection point of the curve in 
Figure 1 gives the best compromise. This corresponds to a period of 
about two seconds for each cycle. It should be noted that a complete 
cycle of the drive mechanism produces two cycles of reactivity 
oscillation. A set of measurements was made on the same sample at 
varying frenuencies of oscillation hut for the same total times. This 
confirmed the cycle time of about two seconds to be proper. 

In order to test the accuracy of the syster , a set of samples 
, was chosen which had varying amounts of boron as an impurity in iron. 
These samples had been measured for impurity both by spectrographic 
methods and by the danger coefficient method in the course of 
investigation rejiorted by Carew and Mow ^ 5 1 • A typical group is 
tabulated in Table 1. 



Table 1 



Samples of Iron Containing Boron 



Seunple 


Weiglit (grams) 


Percentage of Boron 


React ivi 


A 


G.3B1 


0.0001 


200 


B 


3.351 


0.0020 


210 


C 


G.37C 


0.0050 


232 


D 


6.239 


0.C072 


218 


E 


6.261 


0.0250 


267 



i^e) 
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5* Discussion of results* 

(a) Alagnitude of error. 

In order to discuss the error of pile oscillator measurements, 
a sample standard deviation is employed as a basis for determining 
the sensitivity of the technique* The sample standard deviation \rill 
be here defined as the standard deviation of a single integration as 
deteriiiined from several consecutive integrations with the same m 
material sample. The sample standard deviation can be projected to 
determine the sensitivity of a measurement for which the integration 
j>rocess is repeated witi the sajiie material sample* The standard 
deviation of the arithmetic mean of the integrations would diminish 
by the square root of the number of repetitions* A set of typical 
readings for one sample is shown in Table 2 demonstrating the idea* 



Table 2 

Integrator Readings from Integrations on Sample E 
Cycle tiuie: 4*5C seconds; 50 cycles each integration 

Reading (volts) 

29.3 
26*2 
28. C 
29*7 
25*9 

Average = 27.8 

Sample standard deviation = 1*7 or 6% 

Reactivity of sample 270yap approx* 

Magnitude of error 16yU^ approx. 



The error of a measurenent would be diminished by 5 ^ 
approximately if five integrations were made on the same sample. 

The result of numerous measurements such as the one tabulated was 
that the saicple standard deviation for the method corresponded to 
an error of from 20^^ to 10^^. In order to reduce this by a factor 
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of ten vould require deten inin^ the aritlirt-etic uienn from about 
ICC integrations for each sample. If each integration is as rapid as 
two minutes, over three hours would be required for the measurenent 
on a single sample. External influences would undoubtedly come into 
play ir/this time. Employing the danger coefficient technique in a 
repetitive manner over a corresponding elapsed time with the AGN 201 
reactor gave a sensitivi t}»^ of • 

Since the instrumentation syster was capable of accurately 
reproducing the test signal, exoriination of the experimental input 
signal for variations is indicated. The oscillating signal which is 
the input to the integrator was recorded on a strip chart recorder. 
An exact tracing froni a typical run is sho^m in Figure 6. 




Figure 6 

Integrator Input Signal 

The severe variations are apparent. These variations are a 
consequence of the total noise in the system. A discussion of this 
noise is desirable. 

(b) Systeni noise. 

The total noise appearing at the final integrator input will be 
called the system noise. Noise is simply defined as unwanted signal-^8 
The system noise may be traced to three principal sources: reactor 
noise, ion clieuiiber noise, and amplifier noise. 
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Although the equations of reactor kinetics are written for 
continuous functions, the actual physical processes occur in discrete 
steps. At sufficiently high power levels the number of events 
occuring are sufficient to limit deviations from the most probable 
pattern of events. When low powers or small perturbations are 
considered, random fluctuations becon.e significant. 

The reactor noise results from fluctuations in neutron flux 
which are caused by the discrete nature of the neutron creation 
process. It was noted that the delayed neutrons come from radioactive 
decay processes and are therefore subject to fluctuations. Both the 
concentration of delayed neutron precursors and their decay tiroes are 
subject to random variations. It is convenient to express the change 
in flux which is caused by the oscillating sample in the forL of a 
comparison with the variations in flux not caused b}’* the sample 
(i.e. reactor noise). This is simply a signal to noise ratio. 

Frisch and Littler derived an expression for the ratio of 

signal to noise for a single observation of an oscillating sample: 



S(^o) = (Ak)^ 



2TNj. T A^(co) 1 i 



1 + T A^(io) 



where 



=3 change in multiplication produced by the sample 
s 

'T a probability of t!io detector detectinjv a fission 



N 



f 



= mean fission rate 



T 



= total tic.e of the observation 



A(u>) = amplitude res])onse of the reactor. 



Rewriting this as S(co) 




G 



where 



.1 + a“(cj)_ 



X-(CJ) 1 i 



shows the signal to noise ratio with respect to reactor noise 



22 



improves as the square root of the ratio between two j)ower levels. 

Ion chan,ber noise has been demonstrated experimentally “^icj to 
vary inversely as the square root of the ion cheuiber current. Since 
current is proportional to neutron flux, this varies in the same 
manner as the reactor noise. Thus the total noise appearing at the 
linear amplifier varies inversely as the square root of the reactor 
power level. 

The third type of noise is that produced in the electronic 
equipment. This was found to be prin.arily sixty cycle hum. This was 
attenuated in the measuring scheme employed. 

The precision of a measurement using this technique will be 
es-entially detemuned by the system noise described. Since this 
increases as the power level of a reactor is reduced, it will be a 
formidable obstacle to the success of a pile oscillator on a low 
power reactor such as the AGN 2C1. 

In reference an expression was derived which is similar to 

that of Frisch and Littler discussed previously. In this case, 
account was taken of the correlation involved in the type of 
experiment. The expression for signal to noise ratio is 

with F the mean fission rate of the 
reactor, T the total time of the observation on the given sample, 
and factor involving the mean nuc her of neutrons released per 

fission. For a typical measurement (a single integration), this 
expression gives: 

1 C 

F = r.l watts X 3 X 10 f i ssions/watt sec. 

9 / 

= 3 X 10 fissions/sec. 

2 / N 

T = 2.5 X 10 sec, (a typical value) 

5 

I = lO^^k 



S _ / 2F T 

N “ 7T V Xj. 
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Thus for a signal to noise ratio of unity, k rmst be of the 
order of 10 ^ Avhich is approximately ten micro-rho. It would appear 
that the liLiitation on the results has been j)laced largely by the 
reactor noise. 

The effect of the noise on the danger coefficient method is slight* 
In the course of an observation, the noise is involved only when it is 
necessar}’' to determine the reactor period. If the criterion for this 
detennination is the position at which the trace of a linear amplifier 

I 

output on a strip chart recorder crosses a given ordinate, the average 
value of the flux level can be read through tlie noise by drawing in a 
line which is the best estimate of the average signal and accurately 
observing where this line crosses the ordinate. A human operator has 
considerable ability in finding tlie average value of such a signal 
through noise. This technique is illustrated in Figure 7. 




time 

Figure 7 

Time Determination from Strip Chart Recorder Trace 



Thus the statistical advantage of a pile oscillator over danger 
coefficient observations is open to considerable doubt. The noise 
factor may be so disturbing as to outweigh all other considerations. 
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6. Conclusions 



It is concluded that a pile oscillator for the y\GN 2C1 reactor 
utilizing the installed ion chaniber and linear amplifier does not 
give satisfactory results. It is not possible to attain with the 
equipment employed the precision o’;tainable with danger coefficient 
measurements on the some reactor (and using less equipirent) . 

Since the instruiaentation system consistently reproduced a test 
input, output variations with the experimental input are attributed 
to the nature of the information contained in the input signal. The 
reactor noise appears to be the major cause of the variations. 

In order to improve the sensitivity of the pile oscillator, 
either one or both of two things must be done. The quality of the 
input signal could be improved by the use of a better ion chamber 
placed closer to the core. Placement close to the path of oscillation 
might peni.it detection of the local signal to increase signal to 
noise ratio. Another improvement would result from operating the 
reactor at a higher power level than the presently licensed limit. 
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APPENDIX I 



ELEMENT AllY BACKGROUND 

In a nuclenr reactor, the reaction sequence by which neutrons 

are produced and lost is self sustaining if proper conditions exist. 

This is described by the multiplication factor, k, being unity. This 

is the ^'critical” condition for the reactor. If k is greater than 

unity the reactor is supercritical and the neutron population will 

increase. Conversely, if k is less than unity, the reactor is 

subcritical and the neutron population will decrease. 

Physical materials have properties which are important in 

relation to re^^ctor phenomena. The probability that a material 

will interact with specific nuclear particles is described 

quantitatively in terms of the ’’cross section" for the reaction. 

Of particular interest is the probability that an atom of the 

material will absorb a neutron. This is expressed in tern.s of 

the microscopic absorption cross section and denoted by , It 

-24 2 

carries the units of "barns"; one "barn" is 10 cm*". The product of 

the nuiiiber of atoms in a unit volume and the macroscopic cross 

section is termed the macroscopic cross section and denoted by -2T , 

Thus =9 N CF • N is usually expressed in units of atoms per cubic 

a a 

centimeter; therefore has the units of era ^ . The number of 

neutrons in a given volume is teru»ed the neutron density and denoted 
by n. It is sufficiently accurate to state that if the total neutron 
population is in equilibrium, the neutron density (neutrons per cubic 
centimeter) is similarly in eouilibrium. Thus = 0. The partial 

derivative is used since n is also a function of space coordinates. 
However, for a reactor at or near critical, the geor etrical shape of 
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the distribution of n irill not change significantly. 

Since the energy or velocitj^ of a neutron is important in 
describing its behavior, the product nv inhere v is the neutron 
velocity is useful. The term neutron flux is defined as this product 
and is of necessity a function of the energy (or velocity) considered. 
This is denoted by the symbol ^ and carries the units of neutrons 
cm *^sec If the neutrons are in thermal equilibrium with the 
surrounding mediufi:, they are called ^’ttiernal neutrons** and the flux 
is called **thermal flux**. Since thermal neutrons cause most of the 
fissions, knowledge of cross sections at theniial energies is of 
considerable importance. Cross section data are usually given for 
neutrons having a velocity of 22C0 meters per second, the most 
probable velocity per unit velocity for themal neutrons at room 
temperature . 

If a reactor is not critical, there will be a net change in the 
nur.ber of neutrons. This may ))e traced to a col binn^tion of three 
causes: leakage out of the reactor, absorption within the reactor, and 
a gain froi sources of neutrons. The leakage per unit volur.e is 
characterized by a term -DV ^ arising from diffusion theory; D is the 
diffusion constant. The absorption tern; is written ^ utilizing 

the definition of macroscopic absorption cross section. The source 
term may be written sin.ply as S. The resulting equation for the 
neutron balance per unit volume per second becomes 

If . DVV - »« + S (1) 

This is l<no\m as the time-dependent neutron diffusion equation. 

All but a small fraction of neutrons f roi. fission arc given off 
iirniediately with the fissioning of uraniur nuclei. These are tended 
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prompt neutrons. The small remaining fraction (C.0065 of the total) 

arc not emitted immediately but folloif the radioactive decay of 

certain fission product nuclei. A total of six species of fission 

product nuclei iirhich are precursors for the emission of these delayed 

235 

neutrons are known for U .Of the total neutrons emitted per unit 
tice, tlie fraction which is due to each species of delayed neutron 
precursor is denoted ^7 ^ ^ = C.CG65. 

The remaining possible sources of neutrons are from external 

means such as a radium-beryllium source or cosmic rays. 

The total source term may be written: S = S + S, + S for 

p d o 

prompt, delayed, and primary or external neutrons respectively. 

The important characteristics of the delayed neutron species or 
groups are given in Table 3. 

T ab 1 e 3 

Delayed Neutron Characteristics 1 
Group Index (i) Decay Constant, (sec ^ Fractional Yield, 



1 

2 

3 

4 

5 

6 



C.C127 

C.C317 

0.115 

C.311 

1 .40 

3.87 



0.000247 

0.00138 

0.(0122 

0.00265 

0.0(0832 

C.C00169 



Total: 0.00650 ( 



It is important to note that the decay is the inverse of the 
average lifetime and that the actual emission of the neutron can 
occur over a time distributed about the average value. Also the 
fractional yield of a given species is an average and also is 
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subject to fluctuations* 

The discussion of the time behavior of a reactor is greatly 
simplified if the niateri vis in the reactor core are homogeneous. This 
is approxinately the c?.se \fitli the AGN 201; therefore this assumption 
is used. 

The technique of the pile oscillator depends on the measurement 

of the change of the reactor n.ulti])! ication factor resulting from the 

oscillating absorber of neutrons. The effective multiplication factor 

of the re.ictor, k, is defined as the ratio of the number of neutrons 

present in the reactor at a given time to the nunber of neutrons 

present in the ininediate ly previous generation. The concept of generation 

implies a characteristic lifetime for the neutrons which is the average 

time between the “birth” and “death” of a neutron. Since neutrons can 

be lost froi the reactor either by absorption in some nuclei or by 

leakage froii; the reactor boundaries, it is useful to separate the 

effects. Thus a lifetime, can be postulated for thermal neutrons 

in an infinite syster from which tiiere can be no leakage. The 

reciprocal of the macroscopic absorption cross section is equal to the 

mean free path of the neutron; that is, the average distance which the 

neutron travels until it is absorbed. Since 21 has the diii.ensions of 

a 

era \ the mean free path has the proper diiaensions of centimeters. If 

V is the characteristic velocity of a thermal neutron, tlie neutron 

lifetiDie 0 is then . If a factor P. which represents the 

Ko vV^ t * 

probability that a therrial neutron will not leak from tlie reactor is 

introduced, the lifetime for finite geon etry is written ^ P^. 

This concept assumes that all the neutrons diffusing in the reactor 

are thermal neutrons. A more correct expression would include the 
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time a neutron spends slowing down from fission energy to thermal 
energy. Since this time is one or more orders of magnitude smaller 
than the thermal lifetime, its neglect will not be overly serious* 



The factor k is related to the multiplication factor k ^ for 
an infinite reactor of the soi e core material by employing an 
additional factor correcting for the leakage of neutrons during 
the slowing down process. Thus k = k^P^P^. The expressions relating 
these factors to the material and geometry of the reactor are not 
pertinent to this treatment. From the definition of k it is seen that 
if Ic = 1 the reactor is critical; that is, the rate of creation of 
neutrons is equal to the rate of loss of neutrons and the time variation 
of neutron density is zero. The quantity ^k = k - 1 is a useful 
measure of the time behavior of the neutron density. IfAk is not 
zero, the reactor will not be in a steady state condition. The time 



behavior of the reactor may also be related to the term ^’reactivity” 
defined by the relation ^ • A convenient unit is the micro-rho 

vrhere ly^^ = ^ 10 . If k is very close to unity, — 

In a finite reactor, if the neutron density at a given time is n, 
there will be Icn neutrons after the lifetime^ . Thus there will be a 
net increase of -liikn, or a production rate of ^ 

neutrons are proitpt. To correct for the delayed neutrons the factor 

i 

must be subtracted. Finally there will be added those neutrons which 



n, where |3 is the total fraction of neutrons which are delayed. 



come from the decay of the dela^^’ed neutron precursors. If is the 

concentration of a given delayed precursor and X • its 

' 6 

characteristic decay constant, the addition will be ^ 






where the sun. is over the six known delayed neutron groups 
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The e'mation for the neutron balance is thus 



(In ^ k 

It -X" 




( 2 ) 



This equation involves the unlcnowns n and C^. A second equation may 
be written for the fornation of the delayed neutron precursors: 






( 3 ) 



Equations (2) and (3) are the equations of the so called one-group 
space-indej)endent model of reactor kinetics. 

A development sttarting from the diffusion enuation and using a 
more sophist icated analysis which takes some account of the 

slowing down of neutrons from fission energies to thermal energies 
gives the equation: 



[k(l -|3) - l] nv 




^ dt 



P) 



The factor p is the fraction of neutrons which escape absorption 

while slowing down frooi fission to thermal energies. Using the 

definitions: k = 

^ _1 

^ a 

P equation (4) becomes 

t 





where C . = p P^ C . . 

1 ^ f 1 

Thus equation (2^ will be in agreement if the are defined as the 
**ef f ect ive" concentrations of delayed neutron precursors. 
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APPENDIX II 



DESCUIPTION OF THE AGN 201 REACTOR 
^ The AGN 201 is a sirall low power reactor having a licensed 
output of one hundred niiniwatts. With increased shielding this can 
be increased to twenty watts. The general characteristics of the 
reactor are given in Table 4. 



Table 4 



Characteristic s 

Power (licensed) 

^Jaxirau^l Thennal Flux IOC nrw 
Critical Mass 

Core Size (cylinder) 

Moderator (polyethylene) 
Reflector (graphite) 

Shield (lead) 

(water) 

Overall Diameter 
Overall Height 



of the AGN 201 

ICC milliwatts 

6 -2 -1 
4*5 X 10 cm sec 

665 gm (20^ 

enrichment in UO ) 

10^* dia. X 9^^” high 

11 kg 

20 cm thick 
10 cm thick 

55 cm thick (lOOO gals) 
6.5 ft. 

9 ft. 



The ’’glory hole” is a one-inch alun.inum tube which runs 
horizontally through the entire reactor and passes through the 
center of the core. This permits introduction of various materials 
into the region of highest neutron flux. A cross section view of 
the re ictor along a line through the glory hole is shown in Figure 8, 
The variation of then!;al neutron flux with radial position in 
the glory hole is shown in Figure 9. 
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Figure 8 

Reactor Cross Section 
View 
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Figure 9 

Thermal Neutron Flux Through Glorj Hole 
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The rate of absorption of themal neutrons by a material is 
proportional to the then al neutron flux through the cjaterial. The 
themal neutron flux varies with position tiirough the reactor. Since 
k is proportional to the r?.te of absorption of neutrons, k will 
depend on the position of an absorber in the reactor. The relative 
effect of an absorber can be referred to its effect at the center of 
the reactor. This is expressed as an importance function* The 
variation of absorption importance function with position in the 
glory hole is shown in Figure 10. A curve of squared cosine function 
is fitted to the importance function curve over the diacjeter of the 
core. 

1 



36 






1 

c - 

\ 



Absorption Importance Function 




Figure 10 
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